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For a closed bi-partite quantum system partitioned into 
system proper and environment we interprete the microcanon- 
ical and the canonical condition as constraints for the inter- 
action between those two subsystems. In both cases the pos- 
sible pure-state trajectories are confined to certain regions in 
Hilbert space. We show that in a properly defined thermody- 
namical limit almost all states within those accessible regions 
represent states of some maximum local entropy. For the mi- 
crocanonical condition this dominant state still depends on 
the initial state; for the canonical condition it coincides with 
that defined by Jaynes' principle. It is these states which 
thermodynamical systems should generically evolve into. 
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I. INTRODUCTION 

Thermodynamics has been developed at a time, when 
the atomistic nature of matter was not well understood. 
Statistical mechanics tried to fill this gap, when, in turn, 
quantum mechanics as the micro-theory proper was un- 
known yet. Meanwhile quantum theory has envoked 
many changes, but the basic structure of modern theory 
is still close to that formulated by Gibbs some hundred 
years ago and therefore still plagued by old problems. 
Examples are the significance of ergodicity the na- 
ture of irreversibility |^,^, the origin of the second law 
|4-7[ and the relation between physics and information 
|B-10[. In spite of entropy being one of the most fun- 
damental concepts in Gibbsian thermodynamics, there 
have been endless discussions about its nature and its 
precise definition As entropy (-change) is ther- 

modynamically measurable (for any individual system), 
one would like to see its reduction to mechanics still to 
be based on observables. But entropy remains an alien 
concept within this frame: rather than to a point (state) 
it refers to a certain accessible volume in phase space, 
which depends on given macroscopic constraints such as 
energy, volume, etc. Nevertheless, this definition, eventu- 
ally based on ergodicity, seems indispensible to calculate 
thermodynamical state functions. 

On the other hand, one wants entropy to account also 
for the irreversibility, which is absent in the underlying 
microscopic equations of motion. Thus the concept of en- 
tropy appears to be related to the lack of knowledge that 



arises during an evolution too complicated to be traced 
in detail. Some sort of coarse-graining |^ is often en- 
voked in this context, entropy would then be expected to 
increase up to the limits set by those macroscopic con- 
straints. Such an information-theoretical approach also 
underlies Jaynes' principle Jl^ , p^ . This principle aims to 
guide our reasoning in the face of incomplete knowledge 
("unbiased guess"). It is not at all restricted to physics. 
But from this point of view it seems as if entropy and the 
state functions would depend on the physicist's ability to 
minimize his lack of knowledege, which is, of course, un- 
satisfactory. 

Within classical mechanics lack of knowledge may always 
be considered subjective: in principle any observable can 
be known with unlimited precision. This is different in 
quantum mechanics. From the uncertainty principle we 
know that there are always observables that are undeter- 
mined. Nevertheless, in single system scenarios, at least 
one observable can, in principle, be known exactly (pure 
states) , while for compound systems there are even states 
where all observables refering to a specific subsystem are 
unknown, even if some compound observable of the full 
system was exactly predictable (pure but entangled and 
thus locally mixed states [9|JT^). Compound systems 
might evolve from states that contain exact kowledge 
about some observable of each subsystem (pure product 
states) into the above mentioned states, featuring this 
fundamental lack of knowledge . 
So, in the quantum domain we have two possible sources 
of ignorance: one being due to our inability to identify 
the initial state and calculate the evolution exactly, the 
other being intrinsic to the momentary state and thus 
present even for an "infinitely smart demon". 
Here we want to show that in typical thermodynam- 
ical situations the fundamental lack of knowledge by 
far dominates over the subjective lack of knowledge in 
the following sense: Almost all the possible evolutions 
(of which we are typically unable to predict the actual 
one), will eventually lead to states that are character- 
ized by a maximum fundamental lack of knowledge about 
the considered subsystem; this lack is only limited by 
macroscopic constraints. Our considerations will entirely 
be based on quantum-mechanical system-environment- 
scenarios Hiln6|. 
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II. STATE MEASURES AND THEIR AVERAGES 



B. Time and path averages 



A. Definitions 



The entropy we are going to consider is the von Neu- 
mann entropy [ |T4| given by 



S{p) :— — fcTrjplnp} 



(2.1) 



which is invariant with respect to any unitary transfor- 
mation. In a diagonal representation of p this reads 



(2.2) 



where the Wn^s are the eigenvalues of p, which are rou- 
tinely interpreted as the probabilities to find the system 
in the eigenstate \n). The von Neumann entropy is obvi- 
ously related to the Shannon entropy. Another measure 
for the broadness of the probability distribution is the 
purity P 0: 



P(p) :=Tr{p2}. (2.3) 



which is also invariant with respect to unitary transfor- 
mations. In general, these two measures do not uniquely 
map onto each other, but, as will be seen below, under 
some aditional conditions maximum entropy states corre- 
spond to minimum purity states. Because of its less com- 
plicated structure, we mainly consider P and get back to 
S in the end. 

Since we are going to deal with bi-partite systems parti- 
tioned into the considered system, labeled by g ("gas") 
and the surrounding labeled by c ( "container" ) the pure 
state of the full system will be denoted as 



(2.4) 



A,B 



where \A,a){\B,b)) denotes the a-th(6-th) energy eigen- 
state with the energy eigenvalue Ea{Eb) of the gas 
(container) system, with a = 1,2. 



1,2.. .TV: 



B- 



N^{N^) are the respective degrees of degeneracy. 
From any such wavefunction the reduced density oper- 
ator of the subsystem, g, is found by tracing over the 
container system: 



E E ^r(V'.^?)1Aa)M',a'| 



A.A'.B I 



For this state the purity is: 



ABCD abed 



(2.5) 



(2.6) 



Here and in the following a, c. A, C, label the gas, 
b,d,B,D, the container subsystem. (Note that P^ = P'^ 
as the total state is taken to be pure) 



Now, let us consider the time average of some state 
measure M which could be the purity, the entropy, etc. 



A/ :- i / M{\i,{t)))dt 



(2.7) 



Choosing a special parametrization for 1-0), we can con- 
vert the time integral into an integral over the trajectory 
generated by the total system's Hamiltonian for given 
initial conditions. Parametrizing in terms of the real 
{ipi) and imaginary [tp[) parts of its am plitudes in some 
basis {!«)} we can write instead of (2.7) 



M 



J|^(0)) 



(2., 



where \d\ip) \ denotes the "length" of an infinitesimal step 
along the trajectory in Hilbert-space. 
The advantage of this special parametrization derives 
from the fact that the effective velocity 



= ^(^(0)|ij2|V(0)) (2.9) 



is constant on each trajectory and thus independent of 
the actual time t or the special point on the trajectory. 
Hence, the integral (2.8) simplifies to 



L 



10(0)) 



M({v^>a) 



(2.10) 



where L is the length of the path. So, in this parametriza- 
tion of Hilbert space, the time average of any state mea- 
sure equals its path average along any given trajectory. 
We will exploit this relationship to infer the time-average 
(or even the "typical" momentay value) of M for cases, 
in which we know that M is constant over a large part 
of the Hilbert space regio n accessible by the respective 
trajectory, (see Sect. Ill B ). 



III. MICROCANONICAL CONSTRAINTS 

A. Accessible Hilbert space region 

Any thermodynamical system has to allow for a parti- 
tion of its full Hamiltonian in the following way: 

H=:HS + H'' + i with (/)< (i/9 (3.1) 

where and H'^ are the local Hamiltonians of the sys- 
tem and the environment, respectively, with i/^, H'^ — 

0, / is some sort of interaction and the inequality for 
the expectation values has to hold for all states that the 
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system can evolve into under given constraints. If a par- 
tition according to this weak coupling scheme was impos- 
sible, the idea of system proper and surrounding would 
be meaningless. Furthermore, weak coupling is required 
to justify the concepts of extensive /intensive variables, a 
basic feature of thermodynamics (cf. also jl^) 
Making such a partition possible, though, might involve 
a re-organisation of the Hamiltonian. Consider, for ex- 
ample, an ideal gas in a container. in this case is a 
free particle Hamiltonian, / contains the interactions of 
all gas particles with all container particles. In this case, 
the expectation value of / for an energy eigenstate of 
the uncoupled problem -I- H'^ would definitely not be 
small. But this deficiency can be overcome by defining an 
effective local Hamiltonian and an effective interaction: 



HS + Va I' ■= I ~ Va 



(3.2) 



Here, V9 is some potential for the gas particles describ- 
ing the mean effect of all container particles together. 
Substituting the real by the effective parts of the Hamil- 
tonian obviously leaves the full Hamiltonian unchanged, 
but now the partition fits into the above scheme. In gen- 
eral, however, /' cannot be made zero. 
It has often been claimed that a system under so-called 
microcanonical conditions would not interact with its en- 
vironment. This, however, is typically not true (cf. |^,^). 
A thermally isolated gas in a container, e.g., definitely in- 
teracts with the container, otherwise the gas could not 
even have a well defined volume; the only constraint is, 
that this interaction should not exchange energy. If the 
energies contained in the gas and the environment 



are to be conserved, it follows that [p[ 











(3.3) 



(3.4) 



Except for these constraints we need not specify / in more 
detail. All interactions that fulfill this relation will create 
perfectly microcanonical situations. Based on these com- 
mutator relations we find that for any energy eigenspace 
specified by A, B 



Ei^a.^wi' = Eiv^r(o)i 



(3.5) 



is a conserved quantity, set by the initial state. Since we 
want to consider cases here that have zero local entropy 
in the beginning (product states), we get 



ElV'ar(0)|^=ElV'a(0)l'l^?(0)l 



^W^W^B (3.6) 



a,b 



a.b 



where (Wg) are the probabilities of finding the gas- 
(container-) sytem somewhere in the possibly highly de- 
generate subspace characterized by the energy eigenval- 
ues Ej^ This is the only constraint that micro- 
canonical conditions impose on the accessible region of 
Hilbert space. 



B. The "landscape" of in the accessible region 

We are n ot abl e to compute the time average of ac- 
cording to (2.10) for we do not know / in detail, most 
likely we do not even know and L'^ in full detail, 
and even if we did, we could never hope to solve the 
Schrodinger equation for a system involving a macro- 
scopic number of degrees of freedom. This is our subjec- 
tive lack of knowledge. But, as will become clear shortly, 
this lack does not really matter. 
To show this, we proceed as follows: 
The minimum purit y state c onsistent with the micro- 
canonical conditions (3.5, 3.6) and its corresponding pu- 
rity are: 



Y^^\A,a){A,a\ 



min 



(3.7) 



We compute the average of P^ over the total accessible 
Hilbert space region. We will show that this average is 
very close to J^^^n for a large set of cases. Considering 
only these cases, we can conclude that P^ « Pmin 
almost all states within this region. Since the only lo- 
cal state that has P^ = P^j„ and is consistent with the 
microcanonical conditions is p^^j„, all states within the 
accessible region that feature P^ ~ Pmin must yield re- 
duced local states that are very close to pf^m- 
Furtherm ore, as there is no other constraint for the paths 
entering ( 2.1C ) except to stay within the allowed region, 
any path will typically leave the very tiny region of max- 
imum purity where it started if the initial state was a 
product state, and then venture through a region of al- 
most minimum purity along almost all of its length. Thus 
in the long time limit any momentary state will be char- 
acterized by a maximum fundamental lack of knowledge. 
To calculate the Hilbert space average of P^ denoted as < 
pa > we need a parametrization for |?/') — {V'afc^j "0^6^ } 
confined to the allowed region (3.5, |3.6[ ). The Hilbert 
space average can then be written as 



< PS >= 



!PH{M)detTUn' 
J det TY\^d(l),i 



(3.8) 



where (/>„ is the respective set of parameters and T is the 
corresponding fun ction al m atrix. 

According to (3.5) and ( |3.6[ ) the real and imaginary parts 
of the amplitudes that correspond to a degeneracy sub- 
space AB must be parametrized to lie on a hypersphere 
of radius R = \/W^Wg. Thus there is a corresponding 
set of parameters {(jy^^} for each degeneracy subspace 
on which the amplitudes of this and only this subspace 
depend. This means that the functional matrix J-^ has 
block form and its determinant factorizes such that 



< Pf > = 



I P'Uab dot T^^U. 

n^^/det^^^nnrf-^r 



(3.9) 



As a consequence, the average over each term of P^ ac- 
cording to (|2.6|) factorizes and reduces to a product of 
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averages over the degeneracy subspaces that the amph- 
tudes in that very term correspond to. 



C. Hilbert space average of 

It turns out that the averages of the terms of be- 
longing to the same "class" have the same functional de- 
pendence on the degree of degeneracy Nab ■= ^a-^b '^^ 
the subspaces they correspond to and on the probabil- 
ity of those subspaces to be occupied Wab ■= W^Wg. 
Thus, we first have to evaluate the average of each term 
within a given class, add all terms and finally sum over 
all classes. 

There are six classes to be considered: 
I: All four amplitudes correspond to different states 
{{A ^C)V{a^ c)) A {{B ^ D) V {b ^ d)). Then each 
average is a product of averages of products of different 
cartesian coordinates ov er fu ll hyperspheres. As can be 
seen from the appendix ( |7.8| ) those averages vanish. 
II: The "gas indices" (indices refering to the gas subsys- 
tem) of all amplitudes correpond to the same gas state, 
but the container indices to two different container sub- 
spaces {{A = C) A {a = c)) A (B ^ D). They factorize 
into: 



/|^fi'|2det^^^n„rf<^^^^ 



lAB 



,AB 



/det^^^n„rf0. 



AD 



The application of (7 



T 



II 



yields 

_ Wab Wad 
~ Nab Nad 



(3.10) 



(3.11) 



There are N^^NbN'jj terms belonging to that class and 
subspace-combination. Since subspaces that have B — D 
are excluded we find, summing over subspaces: 



ABD 



J2 Nl{N%N't-5BD{N%f)T 



(3.12) 



III: The container indices of all amplitudes correspond 
to the same container state, but the gas indices to two 
different gas subspaces [A ^ C) A {{B = D) A {b ^ d)). 
By repeating the above procedure (II) we get: 



(ml 

Nb 



1 - E(^^ 



(3.13) 



IV: The gas indices of all amplitudes correspond to the 
same gas state and the container indices to the same con- 
tainer subspace but to different container states within 



this subspace {{A = C) A (a = c)) A {{B = D) A {b ^ d)). 
Those terms are of the form 



T 



IV 



JdetT^BYl^^d^AB 



lAB 



(3.14) 



To evaluate this equation one needs (7^) and gets 



IV {W%f{W%YnNAB) _ {WlfiWI^f 



TiNAB + 2) 



Nab (Nab + 1) 



(3.15) 



For each subspace-combination there are Nj^Ng{Ng — 1) 
terms in that class. Thus, summing over subspaces we 
find: 



^iV^BiV^(^^-l)T'^ 



AB 



E 

AB 



(w^ir(w^Br(^B-i) 



Nab + I 



(3.16) 



V: The container indices of all amplitudes correpond to 
the same container state and the gas indices to the same 
gas subspace but to different gas states within this sub- 
space {{A = C) A{a ^ c)) A {{B = D) A {b ^ d)). For 
those terms the above calculation (IV) has to be repeated 
yielding: 



E^" = E 



AB 



{w%Y{w%Y{Ni~i) 

Nab + 1 



(3.17) 



VI: All four amplitudes correspond to the same state. 
These terms read: 



T 



VI 



/i^„i,^|4|det.F-4^n„#r 



/det.F^^n„rf0r 

Using (7J) one gets: 

2{Wi)\WBfnNAB) 



(3.18) 



rpVI _ 



c \2 



T{Nab + 2) 



Nab{N 



AB 



(3.19) 



Since there are NaNb terms in this class for each sub- 
space, one gets: 



Nab + 1 



(3.20) 



If we now finally sum up all the contributions from the 
six classes we get: 



E(^b)= 



Ni 




_j2iwinwi,nNi+Ni,) 



A,B 



N'aNI, + 1 



The Hilbert space average is thus a unique function of 
the invariants W^,Wb, specified by the initial product 
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state, and the degeneracies N^,N^. 
If we, just to get in contact with previous results, ask 
for the average purity of an arbitrary state with no con- 
straints at all we can enlarge the accessible region over 
all Hilbert space by formally taking both subsystems as 
completely degenerate. Doing so, we recover 



< >= 



N3 ^N" 



1 



(3.22) 



as a special case 
If the degeneracy of the occupied energy levels is large 
enough that N^N"^ + 1 



N'^Ng, which should hold true 



for typical thcrmodynamical systems, ( ^.21 ) reduces to 



B 



(3.23) 



The first sum in this expression is obviously exactly P^jj„ 
( |3.7D , so that for systems and initial conditions, in which 
the second sum is very small, the allowed region almost 
entirely consists of states for wich « Pmin- 
second sum will be small if the container system occu- 
pies highly degenerate states, typical for thcrmodynami- 
cal systems, in which the surrounding is much larger than 
the considered sy stem. This is the set of cases mentioned 
already in Sect. [II B. Since the density operator that 
has = P^i„ and = Sf^^^ and is consistent with 
the microcanonical conditions is unique, the density op- 
erators with P^ 

this one and should therefore also have 
more detailed lengthy but straightforward perturbative 
analysis shows that: 



P^j„ should not deviate much from 

Qa A 



(3.24) 



where is a positive function that s cales linearly with 
the system size of the gas system. (3.24) is valid for 
t he th crmodynamical regime, i. e., if the second sum in 
( 3.23 ) is small compared to the first one. 



IV. CANONICAL CONSTRAINTS 

Our approach to the canonical conditions will be based 
on similar technique s as before. The possibility of a par- 
tition according to (3T) is still assumed. But now there 
is no further constraint on the the interaction /, since en- 
ergy is allowed to flow from one subsystem to the other. 
The only constraint for the accessible region therefore 
derives from the initial state of the full system and the 
fact that the probability to find the total system at some 
energy E 



AB i2 
b I 



(4.1) 



A,B/E 



A,B/E a,b 



should be conserved. (Here A,B/E stands for: all A,B 
such that Ea + Eb — E) One could try to repeat the 
above calculation under this constraint, but now it turns 
out that the actual minimum purity is no longer near the 
average purity over the accessible region. Thus, one has 
to proceed in a slightly different way. 
Contrary to the microcanonical case the probability to 
find the gas (container) subsystem at some given energy 
is no longer a constant of motion here. But we are go- 
ing to prove that there is a predominant distribution, 
{^As}' applicable to almost all states within the al- 
lowed region. The subregion formed by these states will 
be called the "dominant region". Once the existence of 
such a dominant region has been established, we can use 
the results from the microcanonical conditions to argue 
that almost all states within this dominant region fea- 
ture the maximum local entropy that is consistent with 
the predominant distribution. 

Therefore, just like in the previous case, our subjective 
lack of knowledge about where to find the system within 
the accessible region should be irrelevant. The reduced 
local state p^{t) as a function of the full state |i/'(i)) 
should always evolve into a state with a fixed probability 
distribution Wj^ and an almost time invariant entropy, 
which is the maximum entropy that is consistent with 
this (canonical) distribution. Nevertheless, the state of 
the ful l sy stem continues to move with the constant ve- 
locity (|2.9| ) in Hilbert space. 



A. Identification of the dominant region 

First, we calculate the size of a region in Hilbert 
space that is associated with a certain energy distribution 
{VFab}. This size will then be maximized with respect to 
the Wab 's under the condition of the energy probability 
distribution of the whole system {W^;} being kept fixed, 
in order to fi nd t he predominant distribution {W^^^}. 
According to (3J) the size of the region asociated with 
the energy distribution {Wab} is: 



A{{Wab}) / det^F^^H' 

AB'' n 



lAB 



(4.2) 



Those integrals are just the surfaces of hyperspheres and 
can be done using the techniques described in the ap- 
pendix: 



AHWab}) ■■^U^AB^^'^'OiNAB) 



(4.3) 



AB 



Here 0{Nab) is the surface area of a 2A'yiB-dimensional 
hypersphere of radius R = 1. 

Instead of maximizing A directly, we choose to maxi- 
mize In A: this is equivalent, since the logarithm is a 
monotonous function. Furthermore we drop all terms 
that do not depend on {W^b} since they are of no rel- 
evance for the maximum and set Nab — 1/2 « Nab, 
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an approximation that is not necessary but simplifies the 
calculation and is definitely valid for large degrees of de- 
generacy. Introducing the Lagrange multipliers {A^}, 
the function we want to maximize with respect to the 
Wab reads: 



Ini - ^ Nab IuWab - ^ ^ Wab - We 

A,B E \A,B/E 

(4.4) 

This maximization is routinely done by solving the fol- 
lowing set of equations: 



ainA 
BWab 



= 



and yields 



W 



Nai 



ab 



E=Ea+Eb 



(4.5) 



(4.6) 



Finally, using (4.1) we find for the Lagrange multipliers 



Ab--^ Ne^ V] Nai 
We ^ 

A,B/E 



(4.7) 



We have thus identified the energy probability distribu- 
tion, which most of the states within the accessible region 
exhibit, i. e., the energy probability distribution of the 
dominant region, {Wj^^}. 



B. Analysis of the size of the dominant region 

So far we have only shown that among the regions with 
given energy probability distribution {W^ab} there is a 
biggest (or smallest) one. But for our argument we need 
to show that this region A'^ is, indeed, extremely larger 
than all the others, that it really fills almost the entire 
accessible region. To examine the size of this region we 
need to know, how the size of a region depends on the 
corresponding distribution {VF^s}, if this distribution 
does not deviate much from the dominant distribution 
{^isl- Therefore we consider Wab ='■ W%g + sab, 
where the e^s's are supposed to be small. For In A we 
then find 



lni = ^ NabHWIb + cab) (4. 

E A,B/E 



with 



A,B/E 



(4.9) 



The latter condition guarantees that Wab — ^ab + ^ab 
still belongs to the accessible region. 



Expanding the logarithm to second order we get: 



In A 



E 

E A,B/E 



^AB 

Wb 



^AB 



W 



AB . 

(4.10) 



Since the expansion is around an extremum the linear 
term should vanish. Indeed, using (4.9) the second sum- 
mation over this term yields: 



E 

A,B/E 



Nai 



<^AB 



W 



AB 



E 

A,B/E 



^EfAB = 



(4.11) 



so that, using (p~q) and (4.7), we finally find 



i'' J|exp- 

AB 



Ha.b/e^ai 
2NabWe 



-^AB 



(4.12) 



i.e. regions. A, that correspond to energy probability dis- 
tributions that deviate from the dominant one are smaller 
than the dominant region, A"^. Since the smallest factor 
that can appear in the exponent of (4.12) for given Ne 



is J the regions A will be extremely smaller already 



for very small deviations if the corresponding A^^'s are 
large. 



C. The canonical distribution 

Finally, to find the marginal, dominant energy proba- 
bility distribution W^ of the gas system one has to sum 
the compound probabilities W^^ over the irrelevant con- 
tainer system to obtain: 



wi = Y^AB = J2 



Nai 
Xe 



^lE 



Ne 



(4.13) 



to 

the state den sities {E^),n'^{E^) and from sums to in- 



Going now from the degrees of degeneracy, N^,iyg, 



tegrals, ( 4.13 ) reads 



W^iES) ^ n3{E3) 



Ea 



W{E) 



n''{E-E3)dE (4.14) 



To better appreciate the meaning of this result we first 
consider the case of the full system energy being exactly 
specifi ed in the beginning {W{E) = 5{E — U)). Then 
( 4.14 ) reduces to: 



W'^{ES 



n'^iU - E9)n3{E9) 



(4.15) 



This is precisely the form postulated in older approaches 
^l| by simply taking the set of local energy eigenstates 
of both subsystems for the full space of the system. In 
this way the authors excluded superpositions or entan- 
gled states, but had to assume perfect ergodicit y, an as- 
sumption that is not needed in our approach. (4.14) is 
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just a combination of such distributions weighted by the 
probabiUty to find the total system at the corrsponding 
energy and then adequately normalized. 
For the remainder we follow those older approaches and 
their assumptions about the spectra of large systems. 
Thus we apply the standard expansion of n''{E — Ea) 
that is supposed to be valid even for large E^, if the 
system is large, 



n^iE -E3)^e 
so that we can write 



lnn'^(£;)- 



r lnn'^(_E)£;9 



(4.16) 



Jea "(^) 

(4.17) 



be 



Assuming the function ^''fj^j'^'^ to 
sharply enough at some value E ~ [/, 



still peaked 
such that 



^lnn'^(£') does not vary much in the range where 
the former function takes on non-negligible values, (and 
again, this range can be fairly wide, if the system is large) 
we get 



W^iEA) oc n^ES)e' 



kT 



(4.18) 



where we have defined -^\nn'^{U) := 
A set-up similar to the one that we discussed to derive 
the canonical distribution, but more spezialzed, has been 
analized by Tasaki |]l6|. His results for which entirely 
different techniques were used and some additional as- 
sumptions had to be made, are in good agreement with 
ours. 



V. DISCUSSION OF THE RESULTS 

A. Attractor states 

The accessible region of the total pure state vector can 
be decomposed into different zones, each of which yield- 
ing local states with a common feature such as purity or 
energy probability distribution. Provided the thermody- 
namic limit (which is defined below) applies to the sys- 
tem, each accessible region is almost entirely filled with 
one dominant zone, corresponding to one special local 
state. This (equilibrium) state can thus be considered 
a local attractor. Hence even for trajectories starting 
in a tiny "non-equilibrium zone" , corresponding to some 
strongly time dependent local state, after some time the 
gas system is most likely to be found in the time indepen- 
dent attractor state, no matter where exactly the total 
system keeps wandering around on its trajectory. 
Those attractor states are as follows: 
(i) Microcanonical conditions: Here we have found 



rmiC. 



Y^^\A,a){AM 



A,a 



(5.1) 



where A labels the respective energy eigenspace, a the 
different eigenstates within an eigenspace, Wa is the 
probability to find the system at the energy Ea and Na 
is the degree of degeneracy of the corresponding energy 
level. 

p^j^ is the state with the highest entropy that the initial 
state can possibly evolve into under strict energy conser- 
vation. 

(ii) Canonical conditions: In this case we have shown 
that 



-HB 



EAae"^IAa)(Aa| 



Tr 



(5.2) 



This state is obviously the same one would have obtained, 
if one had applied Jaynes' principle, taking energy as the 
only relevant observable. 



B. Quantum-thermodynamic limit 

For the above to hold, the system has to fulfill certain 
requirements which then define the (quantum-) thermo- 
dynamic limit. 

For both conditions the weak coupling limit should apply 
to the total composite system, i. e.: 

H^:H3 + H^ + i with (/) < (5.3) 
where and H'^ are the local Hamiltonians of the sys- 



tem and the environment, respectively, with 



0, / is some sort of effective interaction. The expectation 
values are to be taken for energy eigenstates of the local 
uncoupled systems. If the coupling was too strong the 
notion of subsystems would become meaningless, if the 
coupling was too weak, thermalization times may become 
extremely long. 

For a system to fulfill microcanonical conditions there are 
further requirements: 



H\I =0 



(5.4) 



If (5.4) holds, no energy is exchanged between the two 
subsystems regardless of the strength of the interaction. 
Furthermore, the following condition has to be met: 



E 



c\2 



m 



A 



(5.5) 



This holds if the environment system, c, occupies energy 
levels of higher degeneracy than the considered system, 
g, and/or if its energy probability distribution is broader. 
This is likely to be the case, if the environment is much 
larger than the considered system. 

For canonical conditions the full system state density 
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n{E) must be large at those energies that the full sys- 
tem occupies. In addition, the function j^lnn^(£"^) :— 
has to be approximately constant over some energy 
range. Within this range the state density of the envi- 
ronment rf{E'^) has to be higher than the state density 
of the gas system n^{E^) in the corresponding range. All 
this is, again, typical for "large systems" and for the con- 
tainer system being even larger than the gas system. 
In this definition of the thermodynamical limit there is 
no necessity for the systems to consist of many particles 
(cf. [^); nevertheless, the above mentioned criteria are 
most likely fulfilled by such systems. 



C. Ergodicity 

In classical statistics ergodicity (of isolated systems) 
is meant to imply hat the time average (of an individ- 
ual system) equals the ensemble average. Experimentally 
accessible is only the former, but then it becomes ques- 
tionable, under what conditions this equivalence should 
actually hold. And why should the (statistical) entropy 
based on the ensemble distribution have something to say 
about an individual system? 

In the quantum treatment the reduced density operator 
Pg describes, indeed, the individual subsystem in inter- 
action with its individual environment. Obviously the 
underlying quantum nature renders irrelevant questions 
like: How long do we have to wait until the single sys- 
tem has actually visited all those different pure states as 
required by the postulated ensemble-properties. 
Only if one attempts to find a "quasi classical" interpre- 
tation of the reduced von Neumann entropy, the missing 
information about the actual local state may be asribed 
to our "ignorance" 

The fact that an individual system (embedded in a quan- 
tum environment) "simulates" a whole ensemble [p3| may 
be considered a remarkable signature of "quantum par- 
allelism" - a parallelism that has been implicitly antici- 
pated since the beginnings of statistical mechanics. 
This kind of "ergodicity" appears to be a natural ingredi- 
ent of quantum mechanical system-environment scenar- 
ios. While the local entropy is a property of each indi- 
vidual total state vector, for this vector itself ergodicity 
is not needed at all. 



VI. CONCLUDING REMARKS 

In this paper we have considered Hilbert space in quan- 
tum mechanics as the analog to phase space in classical 
mechanics. This analog is supported by the fact that 
in both cases "micro-states" correspond to state space 
points, and their dynamics to deterministic trajectories. 
If the total system is partitioned into two subsystems ( 



the considered or "gas" system, g, and the environment 
or "container" , c) , the reduced local state, p^, for the gas 
subsystem can be calculated from the pure total system 
state l^p). Any such local state has typically non-zero 
entropy and different total system states may very well 
yield the same local system state. This fact, in turn, al- 
lows to define a probabilty distribution of local entropies 
on the space of total pure states. This distribution has 
been the main target of our investigation. 
For composite systems we have specified classes of in- 
teractions, i. e., microcanonical (mic.) and canonical 
(can.) Hamiltonians, respectively, which generate respec- 
tive Hilbert space trajectories. These trajectories can- 
not cross the boundaries of different "accessible regions" 
which are set by the condition type ("mic." or "can."). 
This, again, is reminescent of classical systems for which 
the trajectories may be confined to certain energy shells 
in phase space. 

We have studied the fundamental mechanisms underly- 
ing Jaynes' principle with respect to energy exchange. 
Jaynes' principle, which may be taken as a rule for un- 
biased reasoning, thus coincides with emergent physical 
laws. Similar considerations should hold also for other 
interaction types like particle exchange (grand canoni- 
cal condition). In any case, quantum mechanics should 
play a central role in the foundations of thermodynam- 
ics, not only, e. g., in the low temperature limit. The 
predominance of local equilibrium states does not imply 
that non-equilibrium states were inaccessible or of mi- 
nor interest. In fact any quantum machine, in particular 
the quantum computer, will require significant deviations 
from equilibrium. Machine design must specify means to 
prevent the system from running into the typical equilib- 
rium behaviour. 



VII. APPENDIX 



The integrals that are to be evaluated in Sect. Ill C 
are essentially of the form: 



Z = 



where the u are integers and the z are parametrized by 
the (f> to lie on a hypershere of radius R and dimension d. 
(The zi correspond to the real and imaginary parts of the 
amplitudes, respectively, d is the degree of degeneracy of 
a given subspace and t he pro bability of this subspace 
to be occupied, cf. Sect. Ill B ). The Zg are related to 
cartesian coordinates Xg by: 



R 



(7.2) 



where r is a radial variable. The technique to solve this 
integral, Z, is basically the same as used to calculate 
the surface area of a hypersphere of arbitrary dimension. 
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This surface area is also the special case of Z for ui = 

Um = 0. 

Defining 

Z.:=l e- (-) dr (7.3) 

the product Z2 = ZZ\ may be written as: 

_„2 I Zir\'^i /Zmr\«'" ,r\d-\ 



z. 



m 



R 



{- 

\R 



d-l 



det J- JJ^ d(l)ndr 

n 

(7.4) 



Now, this is an integral over all space, written in angular 
and radial coordinates that can be converted to cartesian 
coordinates yielding: 



Z2 = 



W_dXs 



(7.5) 



Since this integral factorizes completely, both Z2 and Zi 
can be evaluated using standard tables of integrals, and 
depend only on R,d,ui,Um- (Note that Z2 vanishes if 
any of the u is odd.) Z is then found to be 



Z{R,d,Ui,Um) 



Z2{R,d,Ul,Um) 
Zi{R,d,Ul,Um) 



(7.6) 



The averages that are to be computed in Sect. [II C 
of the form: 



are 



A{R, d, ui,Um) 



Z{R,d,uuUm) 
Z{R,d,0,0) 



(7.7) 



They are all invariant with respect to exchange of the 
m's. Here we need only: 

A(R, d, 0, 1) = A(R, d, 1, 1) = 0, A(R, d, 0, 2) = — 

d 

(7.8) 

which could also have been found from simple symmetry 
considerations, and 



i?^r (-) 

A{R, d, 2, 2) = A{R, d, 0, 4) 



4r (f + 2) 



4r (I + 2) 



(7.9) 



All A's are invariant with respect to exchange of the u's. 
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